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We compute directly the entanglement entropy of spatial regions in Chern-Simons gauge 
theories in 2 + 1 dimensions using surgery. We use these results to determine the universal 
topological piece of the entanglement entropy for Abelian and non-Abelian quantum Hall fluids. 



oj 1 Introduction 

The problem of quantum entanglement and its measurement has a long history in quantum me- 
chanics, going back to von Neumann, who introduced the concept of entanglement entropy. The 
quantum mechanical state of a subsystem A is defined by its reduced density matrix obtained 
by tracing out the information contained in £>, the rest of the system. Here A and B are a partition 
of a larger system which is assumed to be in a pure quantum state. The von Neumann entanglement 
entropy Sa (Sb) of region A (B) is defined to be 
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For a quantum mechanical system with a finite (and typically) small number of degrees of freedom, 
the von Neumann entropy is a useful and quantitative way to quantify the entanglement encoded 
in a quantum state. 



The entanglement entropy in a quantum field theory, and for that matter in any quantum mechani- 
cal system with an infinite (macroscopic) number of degrees of freedom, is in general a complicated 
non-local quantity whose properties are not well understood. In quantum field theory, interest in 
the properties of the entanglement entropy arose in the context of finding a possible explanation 
of the Bekenstein-Hawking area law of black hole thermodynamics in terms of quantum informa- 
tion concepts. In a local quantum field theory in d space dimensions, the entanglement entropy 
for a finite region of linear size L scales with the size of the boundary ("area" ) (L/af- 1 of the 
region, but with a non- universal (i.e., dependent on the choice of the ultraviolet cutoff a) prefactor 
without an a priori relation to any general-relativistic quantities. [1] Further studies showed that in 
1 + 1-dimensional conformal quantum field theories this generally non- universal field-theoretic area 
law reduces to a universal In (L/a) size dependence with a universal coefficient equal to c/3, where 
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c is the central charge of the conformal field theory O |3j. Aside from these important results, little 
else is known about the behavior of the entanglement entropy in quantum field theory. 

On the other hand, interest in the behavior of the entanglement entropy in condensed matter 
systems arose in the context of studies of systems near quantum critical points. In that context, the 
concept of entanglement entropy offers a new perspective to characterize the behavior of quantum 
critical points from a unique quantum mechanical perspective. However, although the behavior of 
the entanglement entropy has been studied in a number of interesting quantum critical systems 
[3 HI El El El [8] , its behavior at generic critical points is not yet well understood. Progress on this 
problem is of general interest since deeper understanding of the scaling behavior of the entanglement 
entropy near quantum critical points will shed light on its general structure in quantum field theory, 
and vice versa. In addition, this new way of characterizing quantum phase transitions is of interest 
in the context of current efforts to use such systems for quantum computing. 

It turns out that the theories for which the concept of entanglement entropy is particularly powerful 
are topological quantum field theories. The best understood topological quantum field theory is 
the Chern-Simons gauge theory in 2 + 1 dimensions. [TO] The main purpose of this paper is to 
determine the connection between the entanglement entropy of Chern-Simons gauge theory and its 
topological data. We will work out the entanglement entropy for a general Chern-Simons theory 
with gauge group G and level k on general spatial topologies. We will then apply our results to 
the cases of most physical interest. In addition to their intrinsic interest in topological field theory, 
our results are relevant to the study of topological phases of condensed matter systems whose low 
energy effective field theories are topological quantum field theories. Our results may also have 
relevance in black hole physics in view of Witten's conjecture on the relation between the Bahados- 
Teitelboim-Zanelli black hole [11] of 2+ 1-dimensional gravity and Chern-Simons gauge theory. fT2] 
It might also prove interesting to compare these holographic results with a direct computation of 
the gravitational entanglement entropy using the Chern-Simons description of 2 + 1 dimensional 
gravity. 

Although some of the concepts we discuss here can in principle have wider applicability, in this 
paper we will only concern ourselves with topological field theories in two space dimensions. It has 
been shown fT3| H3] that for a field theory in two space dimensions that is topological in a limit, 
the entanglement entropy for a large simply connected region A of linear size L with a smooth 
boundary (a subset of an effectively infinite simply connected system) has the form 



where a is a non-universal coefficient. This form holds provided the linear size L of the region is 
large compared to any intrinsic length scale of the theory. The universal constant term 7, known 
as the topological entropy, characterizes the topological state and it is a property of a topological 
field theory. For a general topological field theory it is given by [HI E] 



where d{ are the quantum dimensions of the quasiparticles (labelled by i) of the excitation spectrum 
associated with this phase, and V is the effective quantum dimension. [TO] 
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Topological phases in two space dimensions are states of matter which satisfy the following prop- 
erties. They are "liquid" phases, translationally invariant ground states that do not break spon- 
taneously any symmetries of the system. On manifolds with a non-trivial topology (e.g. a torus) 
the ground states exhibit a non-trivial degeneracy which is robust since it cannot be lifted by the 
action of any local perturbation. In these phases the excitation spectrum is gapped. In the limit of 
low energies and long distances the wave functions of a set of excitations, vortices of these fluids, 
exhibit non-local properties that do not depend on the positions of the excitations. The states of 
the excitations of a topological state span a topologically protected Hilbert space, and the rate of 
growth of the dimension of this Hilbert space (as a function of the number of excitations of type 
i) is the quantum dimension d{ of the excitation. In the limit of a vanishing correlation length 
£ — > 0, the effective field theory of a topological phase is a topological field theory. The path 
integral (partition function) of a topological field theory is a topological quantity in the sense that 
it is independent of the metric of the space. The prototype of a topological field theory is Witten's 
Chern-Simons gauge theory of the Jones polynomial [9 1 ITU] . 

The best known and most studied (both experimentally and theoretically) topological phases in 
condensed matter are the fractional quantum Hall (FQH) fluids, incompressible phases of two- 
dimensional electron gases (2DEG) in large magnetic fields (for a review see Ref.[T6l [TT]). Another 
experimentally accessible candidate for a topological phase is the superconducting phase of the 
quasi-two-dimensional strongly correlated oxide Sr2Ru04 which appears to be a p x + ip y supercon- 
ductor and is also a topological state. [HI [19] Ultra-cold bosonic gases in rotating magnetic traps 
have also been conjectured to form bosonic analogues of the fractional quantum Hall states. [20J 

The non-local behavior of the excitations of a topological phase in two space dimensions is closely 
related to the braiding properties of their world lines which, in turn, determine the analytic prop- 
erties of their wave functions. These excitations are generally known as anyons and carry fractional 
statistics. [211 E2] Excitations with Abelian fractional statistics are labeled by one-dimensional 
representations of the braid group and their quantum dimensions di = 1, and their associated 
Hilbert spaces are one dimensional. Excitations with non- Abelian statistics are labeled by multi- 
dimensional representations of the braid group, have quantum dimensions di > 1, and their asso- 
ciated Hilbert spaces are multi-dimensional. Such non- Abelian excitations, and their topologically 
protected Hilbert spaces, are the basis of the concept of topological quantum computation. [231 EH 
[25] 

In spite of its non-local nature, the topological entanglement entropy of FQH states is actually of 
direct physical interest, and important for the characterization of quantum Hall interferometers. 
Quantum interferometers of FQH fluids are devices which in principle can detect the fractional 
statistics of the FQH quasiparticles.[26, 27, 28, 29, 30, 31, 32, 33, 34J It has recently been shown [35, 
[361 137] that a quantum point contact in a non- Abelian quantum Hall state actually can act as a 
quantum disentangler. This is possible because point contacts in a quantum Hall system essentially 
consist of places where the tunneling matrix elements between the edge states is non-vanishing. 
From the point of view of the bulk FQH state this is a non-local connection which disrupts its 
quantum correlations. Remarkably, Fendley and coworkers [35 1 136 1 [37] found that the change of the 
entanglement entropy induced by the point contact is equal to the change in the Affleck-Ludwig 
entropy [38] of the edge states of the FQH fluid. Thus, the topological entropy is a quantity of 
interest in the present effort to develop interferometers for quantum Hall quasiparticles. 
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In this paper we investigate the universal properties of the entanglement entropy of Chern-Simons 
gauge theory for a general compact gauge group at an arbitrary level k. We give explicit re- 
sults for SU(2)k and for several coset conformal filed theories of interest for applications. We use 
our results to compute the entanglement entropy of fractional quantum Hall topological fluids by 
computing this quantity directly at the level of the effective topological field theory. By working 
directly in the topological limit we obtain directly the O(l) term in the entropy, the topological 
entropy. All other size-dependent terms, including the "area term", become zero in this limit. 
Naturally, size-dependent terms will arise if (irrelevant) corrections to the topological action, such 
as Maxwell/ Yang- Mills type terms, were to be included. Throughout this paper we will use the 
path integral representation of Chern-Simons theory. To this end we will adapt the seminal results 
of Witt en [9} [TO] for the partition functions of Chern-Simons gauge theories to the computation of 
the topological entanglement entropy. We use the standard "replica" approach to compute the 
entropy [321 13] • This requires to understand what is the 3-manifold resulting from gluing n copies 
of the system in a suitable fashion [3], needed to compute the entropy for a number of cases of 
interest. The key aspect of our approach is the identification of a suitable configuration of Wilson 
loops for each case of interest and to compute it by reducing it to already known cases by using 
surgeries. Alternatively, it is also possible to use a more conventional approach using the wave 
function of the Chern-Simons gauge theory [10J. This approach is technically more involved and 
will only be discussed briefly. 

We consider first the case of a surface of genus zero (a sphere or a disk). 1 For a simply connected 
region we compute the topological entropy of the vacuum (ground) state of the Chern-Simons theory 
on a sphere, and recover the result obtained by Kitaev and Preskill[13j, and of Levin and Wen|14|, 
i.e., eqs. @ and @. Next we generalize these results to the case of manifolds with non- vanishing; 
genus (mainly a torus), which have a finite-dimensional topologically protected degenerate vacuum 
sector. Here we also consider the entanglement of multiply connected regions. We find that the 
entanglement entropy of a simply connected region is independent of the genus of the manifold, 
even if the vacuum sector is degenerate. In the case of a multiply-connected region, we find that the 
entropy scales linearly with the number of components of the observed region only if the vacuum 
sector is non-degenerate (genus zero). However, if the manifold has a non- vanishing genus, and thus 
has a degenerate vacuum sector, in general the entanglement entropy of multiply-connected regions 
is different for different states in the vacuum Hilbert space. In other words, the entanglement 
entropy, aside from the purely topological entropy, has additional contributions that depend on the 
choice of state i.e., on the coefficients of its wave function and on the representation carried by 
the state. We also compute the entanglement entropy of a simply connected region with several 
quasiparticles, i.e. operators represented by punctures carrying non-trivial representations. In this 
case we find that the entanglement entropy generally depends on the conformal blocks in which these 
operators can fuse, and hence depend explicitly on the structure of the fusion rules. These results 
indicate that measurements of the entanglement entropy can, in principle, be used to determine 
the full structure of the underlying effective topological theory. Finally, we apply these results to 
the computation of the entanglement entropies of fractional quantum Hall fluids. Here we derive 
the modular 5-matrices for several coset CFTs needed to compute the entanglement entropies for 
non-Abelian FQH states. 

1 As this manuscript was being completed we became aware of the very recent work of K. Hikami [40] . who 
calculated entanglement entropies in SU(2)k theories on the sphere using a skein relation approach. Our results 
agree with Hikami 's where they overlap. 
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This paper is organized as follows. In Section[2]we set up the calculation of the topological entangle- 
ment entropy 7 in Chern-Simons gauge theories. In Subsection |2.1| we show how the computation of 
the entropy can be carried out using the methods developed by Witten [9] and use them to compute 
the entropy for the simplest cases, a simply connected region on a sphere (Subsection |2.2|) and a 
torus (Subsection |2.3| ), and multiply-connected regions on a sphere (Subsection |2.4[ ) and on a torus 
(Subsection |2.5| ). In Section [3] we present a calculation of the entanglement entropy in the presence 
of punctures (i.e. quasiparticles) carryi ng d ifferent representation labels. Here we discuss the case 
of four quasiparticles on S 2 (Subsection 3.1) and discuss two different cases, paired and not paired, 
as well as three quasiparticles on S 2 (Subsection 3.2). In Sect ion [I] we present the calculation of 
the entanglement entropy for both Abelian U(l)k (Subsecti on |4.l| ) and non-Abelian FQH fluids, 
in terms of coset Chern-Simons gauge theories (Subsections |4.2| and 4.3). Section [5] is devoted to 
the conclusions. The hydrodynamic, Chern-Simons, description of the FQH fluids (both Abelian 
and non-Abelian) is summarized in Appendix A, and the calculation of the modular 5-matrix in 
Appendix B. 



2 Entanglement entropy and Chern-Simons gauge theory 



In this paper, we will consider entanglement entropy in Chern-Simons theory in three dimensions. 
As we will argue in what follows, the entanglement entropy may be obtained by computing a Chern- 
Simons path integral on certain 3-geometries, which we systematically obtain through a 'gluing' 
procedure. To see what this procedure should be, we review here the conceptually simpler case of 
a scalar field theory. The 2-dimensional case was described by Calabrese and Cardy.[3j Consider 
a spatial domain which we slice into two, labeled A and B. These regions may be connected or 
not, simply connected or not. We label the interface between A and B by / = dA — dB, which 
may in general consist of several components. We label the degrees of freedom as <j). To make 
the discussion more straightforward, we will consider the corresponding finite temperature density 
matrix 

= j n r )i e ~ SE n 5 °) - )] 5 [<p(x, p) - WW] , (5) 

X,T X 

where we specify the state by a spatial configuration at r — 0,/3. In this language, a trace is 
obtained by path integration over (fo and (pp. Having split the spatial domain into pieces, we may 
then obtain the reduced density matrix pA by tracing over £>, 



PA [{<po(x)},{<pp(x)}\x e A] = J y\[ [d<po(£)d<p^ 

(6) 



KxEB 

The entanglement entropy will be obtained by a replica trick, 



S A -tip A In p A -^-tipA (7) 
an n=i 
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(we expect that tr//^ will have a unique analytic continuation in n for n > 1). Finally, trp^ is 
obtained by taking n copies of pa and 'gluing' them together appropriately 

= I n {n^o^)^?^)] n * - n 5 K^) - ^ 

This path integral may be interpreted as a scalar field theory defined on a glued manifold, of the 
form displayed in Figures |1|2| 



V 




Figure 1: Conceptual picture of pa- The 
trace over B corresponds to gluing r = to 
r = j3 in the B region, leaving a cut open in 
the A region. 



Figure 2: trp\ is obtained by gluing three 
copies of the diagram in Fig. [I] back to back 
along the cut in the A region. 



Now, in fact we are not really interested in the entanglement entropy obtained from the finite 
temperature density matrix. Instead, we would like to pick a pure state (of the whole system); this 
may be achieved here by taking (3 — > oc. In this limit, the system will project down to the ground 
state. There is a subtlety here, that will in fact arise in the Chern-Simons theory (or generically in 
any topological field theory), in that the ground state need not be unique. Thus, the procedure we 
have outlined is not powerful enough to select a particular degenerate pure state. In the particular 
case of Chern-Simons theory, we will take the above construction as indicative that we should 
consider the Chern-Simons path integral on the glued geometry; in this construction, it is clear 
how to make the choice of a pure state, as we will detail a little later. 

In the case of Chern-Simons theory, we can formally perform the above construction by identifying 
the Chern-Simons wave functional; assuming holomorphic factorization, this may be written as a 
WZW path integral 2 

{Bi\{Ai\*) ~ J [dgA,idgB,i]e~ klA ^ 9A,i ^~ klB ^ gB, ^~^ ^ A tT ^ i9 ^ id9A,i ~^ tr ^9B^i d 9B,i ^ (9) 

Here, 1(g) is a WZW action. The expression Q should be interpreted as a sum over histories with 
a spatial section of fixed topology. Here, we have split the integral into contributions of fields in 
regions A and B. Formally, trp^ may then be constructed by gluing together suitable such factors, 

/n 
]J[dii{A k )di^(B k )}(B 1 \{A 1 \^){^\B 1 }\A2}(B2\{A 2 \^)(^\B2)\A 3 ) . . . (S n |(i n |*)(*|S n )|^i). 

k=l 

(10) 

2 Written in this form, the gauge field measure includes a factor / tr ^ A . 
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This can be interpreted as Chern-Simons theory on a glued 3-geometry. This 3-geometry will be 
determined by a choice of spatial topology (that is, a Riemann surface S of genus g) and a choice 
of cutting into A and B regions. The original 3-geometry (before gluing) may be mapped to a solid 
geometry E, consisting of E and its interior. It is well known that the Hilbert space Ti^ of this 
theory is accounted for by the appropriate conformal blocks of the corresponding WZW conformal 
field theory. For example, for the sphere S 2 there is a unique state, while for the torus T 2 , the 
various degenerate states may be obtained by placing Wilson lines in representation R along the 
centre of the solid torus. For higher genus, we can consider the various conformal blocks directly. 
Thus, the choice of pure state is made here by a choice of conformal block (equivalently, for the 
torus, a choice of Wilson loop). 



2.1 Modular Properties 



Indeed, it is well known that the states of a Chern-Simons theory are accounted for by the con- 
formal blocks of a conformal field theory. As a result, the Chern-Simons states may be identified 
with characters. The modular 5-matrix of the conformal field theory will then enter in calculations 
of Wilson loop observables in the Chern-Simons theory, as was exemplified by Witten[9]. Con- 
sequently, the entanglement entropy will generically depend on matrix elements of the modular 
5-matrix. Given a set of characters %(r) of a CFT, one writes 

X (-1/t) = S X (t), (11) 

which should be understood as matrix multiplication. The characters are indexed by a set of 
quantum numbers, which in the case of afflne algebras can be taken to be representations. In this 
paper, we will not need to specify the precise Chern-Simons theory (that is the calculations are 
valid in general), although physical applications will imply a choice. An example of interest is 
SU(2) k WZW, in which representations Rj are labelled by a half-integer j = 0, 1/2, fc/2, and 

xf\-lM = J2 sik) /xf(r), (12) 
f 

where 



° J 



k + 2 



sin 



7T(2j + l)(2/ + l) 

k + 2 



(13) 



More generally, the 5-matrix is assumed to be unitary. Thus, we have 

<^(<St)/ = o?, (14) 
while applying S twice corresponds to charge conjugation 3 

(<S 2 )/ = Ci j = Si. (15) 

Also of importance are the fusion rules for two representations Ri x Rj. The multiplicity of rep- 
resentation Rk in the fusion is denoted by Nij k which are related to the modular 5-matrix by the 
Verlinde formula 

N ^ k = ^ s/s/(s-V (lg) 

<->0 



3 The notation j refers to the conjugate representation to that labelled by j. 
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The quantum dimension is defined as 



For SU(2) k , the quantum dimensions are 



dj = % (17) 
<->o 



sin 



7T(2j + l)\ 

k+2 J 



dj = ^ ^. (18) 



In Section [I] we give generalizations of these formulae to other CFTs of interest. 
We note that the unitarity condition implies 



(So )"^ 4 /Z>il 2 = 2 >- (19) 



In the Chern-Simons theory, we will be led to evaluate the partition function on various 3-geometries 
with Wilson loops. These can be systematically computed by a series of "surgery" operations [9j. 
The result of these computations is that the partition functions depend on various matrix elements 
of modular matrices. For example, 

Z(S 3 ,Rj)=S i, (20) 

where the notation on the left means the Chern-Simons partition function on S 3 with a Wilson 
loop in representation Rj (where j is an index labeling representations). 

Another basic result that we will use repeatedly applies to a 3-manifold M which is the connected 
sum of two 3-manifolds Mi and M 2 joined along an S 2 . We have (Eq.(4.1) in [9]) 

Z(M) • Z(S 3 ) = Z(Mi) • Z(M 2 ). (21) 

This result relies crucially on the fact that the Hilbert space for S 2 is one dimensional. Similarly, 
using the same reasoning, we can deduce that if M is Mi and M 2 joined along n S^'s, 



Z(Mi) • Z{M 2 



m) - ( 22 ) 



To demonstrate this, we note that the path integral on a connected sum of Mi and M 2 (that is, Mi 
and M 2 joined through £ 2 's) can be thought of as the overlap of states |xi) and \x2) eacn defined on 
the interface S 2 , Z[M\ + M 2 ) = (xi|X2)- Because the Hilbert spaces involved are one- dimensional, 
we can insert a state in between which topologically corresponds to capping off the connection (that 
is, we sew in half of a 3-ball onto Mi and half of a 3-ball (with opposite orientation) onto M 2 ). 
Thus, 

7(M M M x (xilX3)(X3|X2) Z(Mi) • Z(M 2 ) 

z(Ml + M2) = (xsixs) = zm ' (23) 

the latter equality obtaining because (X3IX3) corresponds to a 3-sphere. We can repeat this con- 



struction, capping off each join, to obtain the more general result Eq.(22). In the following, this is 
the general surgery operation that we employ repeatedly. 
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2.2 S 2 with one A-B interface 



Let us begin with the simplest case, in which the spatial topology is a 2-sphere. The Hilbert space 
on S 2 is one dimensional, and so there is only one choice of state. The 3-geometry is the 3-ball 
shown in Fig. |3J If we take the A and B regions to be connected, then they are disks. To 

b 





Figure 3: Shading implies a solid 3-ball. With a one-component interface, the A and B regions are 
disks. It is useful in the following constructions to view the 3-ball as a disk rotated about an axis 
passing through the origin, as shown at right. 

construct trp^, we glue 2n such pieces together. In Fig. |1J we show how to systematically perform 
this gluing. We have drawn the n — 2 case explicitly, but it is not hard to generalize to higher n. 
In the figure, we have used 1 and 2 to label \^)\ and 3 and 4 for \ijj)2 and (V>|2- The four 




Figure 4: For spatial topology S 2 with one interface component, we explicitly show the construction 
of trp^. The overall manifold is generated by four pieces of disks glued together one after another 
and rotated along the same axis as in Fig. [3| 

slices form four 3-balls (or as shown, rotated disks); when glued together to form trp^, we find an 
S 2 , rotated about the axis, which has the topology S 3 . One can easily check that for higher n, we 
obtain the same result, the S 2 being obtained by sequentially gluing 2n disks. Thus we have the 
normalized trace, 

tr />V,i) _ Z(S S ) 



= (Z(5 3 )) 1 -" = (So 



0\l-n 



(24) 



where we have used formula ( |20| ) given above. Finally, using eq. Q, we obtain 

sf'^lnSo . 

This result applies to any Chern-Simons theory. Since <So° — we recover the known result [E 
Q3] for the topological entropy for a simply connected region of a sphere S 2 (or a disk) in terms of 



(25) 
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the effective quantum dimension V: 



In the case of U(l) m , we have 



and hence 



For SU(2) k , we obtain 



S { f' 1} =fn<S ° = -InV. (26) 



«5o° = . = -j=, (27) 



S (SM) = _ ln ^_ (28) 



For fc = 1, 2, 3, this evaluates to 

S A = ~lnV2, -In 2, -ln^/>/5 + 5^ , (30) 
respectively. For applications of these formulae to physical models, see Section [4j 

2.3 T 2 with one component A-B interface 

The Hilbert space on T 2 is isomorphic to the space of integrable representations Rj of the Kac- 
Moody algebra. These states are generated by doing the path integral on a solid torus with Wilson 
loop in representation Rj lying along the non-contractible loop at the centre. We will consider first 
a slicing of the torus into A and B regions such that there is a single connected interface component. 
To define the entanglement entropy, we must choose a pure state, and here we have a choice. To 
begin, let us first choose the trivial representation (equivalent to no Wilson loop). It is useful to 
consider the solid torus as a solid ball with a handle attached in the A region, since we have already 
studied the solid ball in the previous subsection, and to compute trp^ here, we need to follow that 
analysis and also keep track of the gluing of the extra toroidal fixtures. Note that a solid torus 
can be thought of as D2 x 5 1 , and two copies glued together (with opposite orientations) gives an 
S 2 x S 1 . The result of the gluing for n — 2 is shown in Fig. [HJ Thus the resulting manifold will 
be the connected sum of an S 3 and n S 2 x S^'s joined along n £ 2 's. Thus, applying Eq.(22) for 
M\ — S 3 and M2 n disjoint copies of S 2 x S 1 , we obtain 

ggW) _ 1 Z(S S )Z(S 2 X S'T _ (7( o^l-n _ re 0,1-n „ n 

(tr^d* i>) n " Z(S* x zWf -W S » ■ ^ 

The first factor after the first equal sign comes from the normalizing factor (for n = 1, the topology 
is just S 2 x S 1 ). We note that this result coincides with the S 2 result. As we shall see, the 
commonality of these two examples is that the interface is the same; the topology of the A and B 
regions themselves does not contribute. 

It is simple to repeat this construction for other pure states, that is including a Wilson loop in 
representation Rj inside the solid torus. This Wilson loop is as shown in Fig. ??. In the gluing 
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S 3 



Figure 5: The space X n= 2, obtained by gluing n = 2 copies of the space shown in Fig. ??. The 
result is an S 3 joined to two copies of S 2 x S 1 along S 2 's. For general n, the glued geometry X n 
consists of an S 3 joined in this way to n S 2 x S^'s. 



above, we will now have a D2 x with Wilson loop in representation Rj glued to a D2 x of 

opposite orientation with Wilson loop in representation Rj (the conjugate state), as indicated in 
Fig. [5j Thus, we have 

tip" 

>0 



tip 



A(T 2 ,Rj) 



Z(S : 



3\n 



= Z{S 3 ) l ~ n = (So ) 1 '". (32) 



Z(S 2 x S 1 1 R jl R j ) n 

In fact this result can be generalized further, to any pure state = J2j ^jlA?)' 

tr /V 2 ,</0 _ ^„ ,/,.... ^^^^3 • • • i>iJ>* h Z{X n , R h ,t h ,...) 



( tr PA(T 2 ,V)) n 



Eil^-I^xsi,^,^))" 



(33) 



where we have denoted the glued 3-geometry as X n \ the Wilson loops Rj and i?j are located along 
the toroidal fixture. Performing surgeries as in Eq.([22|) gives 



tip 



A(T 2 ,ijj) 



(trp A ( T 2^) 7 



E,,, 2 .... M^ Ml ' • • M*n Uk Z (S 2 x S \ R k , R k ) Z(S 3 )^ 
Y J M 2 Z{S 2 xS\R J1 t 



(So 



0\l-n 



(34) 

So we conclude that, at least for this slicing into A and B regions, the entanglement entropy is 
insensitive to which degenerate pure state we consider. This statement is not generally true, as we 
will see. 



2.4 S 2 with two- component AB interface 



Next let's study the case where A and B meet at an interface with two components. For S 2 , the 
only distinct choice is to have two disconnected B regions. If we think of this as two solid 3-balls 
joined together, then we can do the gluing for each 3-ball separately, and then account for the 
joining. The result is a pair of S^'s, joined along n — 2 S^'s, as indicated in Fig. [6] for the n — 2 
case. Thus we find 

_ (Z(S*))-"Z(S S )Z _ m<?3 »2(l-») - K <M(l-n) ft* 
71 \n - C7(Q3\\n ~~ V Z W ) — 1^0 J '■ (ot>J 

(trp A (52,2)) {Z{S 6 )) n 
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Figure 6: S 2 with two-component interface. For n = 2, gluing two copies of Fig. ?? together gives 
a topology that can be thought of as two 5 3 's that are joined along n — 2 S* 2 's. 

It is not difficult to envision the generalization of this result to an M-component interface on S 2 , 
and we find 

= MlnS . (36) 



2.5 T 2 with two- component AB interface 

In the case of T 2 with a 2-component interface, there are a number of new choices to be made. 
There are essentially two distinct ways to slice the spatial surface, which we consider in turn. 

2.5.1 T 2 : Disconnected B regions 



The first possibility is shown in Fig. [7| Here, since we have learned that it is the features of the 





Figure 7: The first of two ways of slicing the toroidal space with a two-component interface. B has 
two components. This may be thought of as two 3-balls joined by two tubes. 

interface that matter to entanglement entropy, we expect that this would give the same result as 
in the last subsection. The glued geometry is shown in Fig. [8] for n — 2. Note that the Wilson 
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Figure 8: The glued geometry for n — 2 in the case of two disconnected B regions on a spatial 
torus. 



loops are located as shown in the figure. For this case, we find 



tTp A(Ti,2,R) 1 ((Z(S 3 )fZ(S 2 xS\R,R)r (7(q W(l-n) / c O^l-n) 

(37) 

which indeed is the same result, for any representation i?, as for the spherical topology with two 
interface components. 



2.5.2 T 2 : Connected B region 



The second possibility, shown in Fig. [9j will present new complications. The new feature here is 
that the Wilson loops thread through the interface between the A and B regions. We will find 
that this leads to a dependence on the representation in the entanglement entropy. Again, this 




Figure 9: The second of two ways of slicing the toroidal space with a two-component interface. B 
has a single component. Again, this may be thought of as two 3-balls joined by two tubes, but a 
Wilson loop threads the interface. 

computation can be thought of as a pair of 3-balls connected by tubes. Upon performing the gluing, 
we obtain a pair of S^'s, connected along 2n S^'s, with Wilson loops routed as shown in Fig. 10 
We find 

tr K(T 2 ,2^)) = Z{S^R 3 )-^Z{S^R 3 )^Z{S^R 3 f = z{s3 k )2{1 _ n) = (SQj) 2(i-n) (3g) 
tr (PA(r2,2,^)) n Z(S 2 xS\Rj,tj) n ' J 
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Figure 10: The glued geometry for n = 2 in the case of a single connected B region on a spatial 
torus. 



To obtain this result, we have generalized the formula (22) to the case where we also have a Wilson 
line. We can do this because the Hilbert space of Chern-Simons theory on S 2 with two complex 
conjugate charges is one dimensional, just as S 2 with no punctures. As we explained earlier, the 
one-dimensionality allows us to cut and glue in half-S' 3 caps; the same formula works for the half-S' 3 
with two conjugate punctures on the surface of the cap. The connection of these punctures inside 
the S 3 forms a Wilson loop. 



Thus, in the numerator of Eq.(38), the first factor corrects for the inclusion of the endcaps (a total 
of 4 half-5 3 's per tube), the second factor comes from each of the 2n capped tubes (which are just 
£ 3 's with a Wilson loop Rj) while the last factor comes from each of the two "large" S^'s that each 
have been capped 2n times. The routing of the Wilson loops through the original 2n tubes is such 
that after all of this surgery, there is a single Wilson loop in Rj on each of these large S 3 's. 

It is then straightforward to show that for any state = J2j 

^Pa(^ M ) = KA(ZijMiKB\RM&\) n 
( tr PA(T 2 ,2,V0) ( tr PA(T 2 ,2,V0) 



( tr PA(T 2 ,2,V>)) n 
_ Ylp=l J2j p \ rt l ) 3p\ 2tr (pA(T 2 ,2,R jl ) • ' ' ^4(T 2 ,2,£ Jn )) 

( tr PA(T 2 ,2,V>)) n 
( tr PA(T 2 ,2,V0) n 

(E.-IV'il 2 )" • ( j 

Here we have used the fact that each S 2 we cut along should have total charge zero, otherwise the 
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path integral vanishes. We thus obtain the entropy 

S A (T^2,i>) = E [2|V>i| 2 ln<V - l^| 2 ln|Vi| 2 ] • (40) 

3 

Since we can interpret |^j| 2 as a probability pj, the second term has the familiar form —pin p. More 
precisely, note that this can be rewritten 

j 



In- 



<3 



(41) 



These calculations can be generalized to higher genus spatial surfaces, using similar techniques as 
we have displayed here. The entanglement entropy is sensitive to the topology only in cases where 
we choose carefully the interface between the A and B regions. Finally, we note that Eq.(41) is 
indicative of a more general result that says that the entanglement entropy depends on the number 
of interfaces, the states and how they fuse, and their quantum dimensions. Notice that in Eq.(41) 
the quantum dimension dj appears squared. This is so because there are two interface components 
in this case. In general there will be a factor of a quantum dimension for each interface component. 
However, in general the entanglement entropy will also depend on the non-universal amplitudes 
in which the state appears in the conformal block. Hence in general the entanglement entropy 
depends on the universal properties of the topological field theory and on the specific form of the 
state. 



3 Quasiparticle Punctures 



It is also of interest to consider entanglement entropy in the presence of quasi-particles. These will 
correspond to punctures on the spatial surface, and to each puncture we associate a representation 
Rj. Here we will consider just the simplest possibilities, but in so doing, we will explore how to 
write the entanglement entropy in a more convenient basis, that of conformal blocks. 




Figure 11: The wavefunctional with quasiparticles is related to a history with Wilson lines coming 
in from t = — oc to punctures on the spatial surface. On the left, time runs vertically to a spatial 
slice. On the right is the equivalent 'radial time' view. 
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3.1 S 2 with four quasi-particles 



A simple case is S 2 with four quasiparticles. We will focus on SU(N) k , where N > 2 and k > 2, 
with punctures carrying 2 fundamental and 2 antifundamental representations on S 2 . We will use 
a and a* to denote fundamental and antifundamental representations, respectively. We may think 
of these punctures as being connected by oriented Wilson lines that extend into the interior; in this 
sense, they correspond to timelike Wilson lines extending in from t = — oo. We note that these 
lines may braid. So we should expect that entanglement entropy may sense this braiding. Now let 
us consider the entanglement entropy, where we simply divide the sphere into two halves. There 
are actually several distinct cases to consider. If the interface between A and B cuts one Wilson 
line (that is, A contains one puncture, say a, and B contains three, the result will be 

S A = In So" (42) 

If A and B each contain two punctures, there are two possibilities. The first possibility has a and 
A* punctures in both A and B regions; in this case, the Wilson lines could connect 6ia®*a and 
aB®*B, or they could connect and a#c¥^. The second possibility is that A contains two cVs 

and B two A*'s; in this case, there are two possible connections of Wilson lines, and these differ by 
braiding. We will see that these choices correspond to choices of conformal blocks in the fusion of 
a with A*. 



3.1.1 B with a and a* 



For k > 2, the Hilbert space on S 2 with 2 pairs of a and a*'s is two dimensional. Let's pick the 
two linearly independent states as follows. After gluing along B, 4 there are four types of density 





i* 2 > 



Figure 12: The two states represented by Wilson lines connecting punctures. 



matrix obtained; these are shown in Fig. 13, If we begin with a pure state \<j>) — a\<j>\) + & 1 02 ) 5 
we have pa — aa*pn + a&*pi2 + Q*bp2\ + &6*/?22- Gluing n copies of these together to form trpA n 
gives rise to an S 3 made from all the possible combinations of pij's. To compute trpA n , we need to 
identify the Wilson loops formed in each case. Each combination contains a number of fundamental 
Wilson loops, each of which contributes a factor ^P-q = d&. For each appearance of pn, there will 
be two such factors. For each factor of P12 or p2\ there will be one such factor. Finally, factors of 
P22 do not increase the number of loops; however, when the contribution to t?pA n is n factors of 



The gluing map is taken to identify punctures, and thus connects the Wilson lines. 
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P22>> there are two loops. Thus, we arrive at 



o o 



E 

j,k,l 



j\k\l\{n - j - k - l)\ 



(aa*) j (ab*) k (a*b) l (bb*) n ~ j ~ k ~ l d 2 J +k+l + (bb*) n (d% - 1) 



= [aa*d| + ( ab * + a*b)d & + bb*] n + (bb*) n {d% - 1). 
After normalization, 



(43) 



(44) 



/ c 0\1 — n 

Thus we find 
where 



aa*d| + (a6* + a*b)d a + 66* 
_(aa* + 66*)d| + (a6* + a*6)d&_ 



+ 



66* 



_(aa* + 66*)d| + (a6* + a*6)d&_ 



[4 - 1] 



S A = ln5 -AilnAi-(4-l)A 2 lnA 2 , 

|ado, + b\ 2 



(45) 



Ai 

A 2 



|ad a + 6| 2 + (d| -1)|6| 2 ' 

\ ad& + b \2 + {d 2 & _ im 2- 



(46) 



We notice that while the definition of \<pi) and | ^2) makes the calculation transparent, they are not 
orthonormal. In fact, we have 



da 1 

1 (lev 



If we define 



d\-i 



(47) 



(48) 
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we can show that the new states are orthonormal. In fact, the new states ](/)[} and l^) correspond 
to the conformal blocks associated with the trivial and adjoint representation respectively, which 
appear in a x A*. We have just calculated the fusion matrix. |<^_) and Ify) are the conformal blocks 
in one channel, while similarly — i_|0 2 ) and / — = =(— 102) should be the blocks 



in the other channel, using their relation, we can easily get the fusion matrix 



1 [ a ex. J 



^6: 



1 



1 



1 



(49) 



where, is (as before) the quantum dimension of the fundamental representation a, and \ j d| — 1 
dn is the quantum dimension of the adjoint representation 6. 



In the conformal block basis the amplitudes become 




(50) 



In terms of the wavefunction in the orthonormal conformal block basis, 



,/|2 



U\2 



Ai = 



z'| 2 + l^'l 2 ' 



A 2 = 



(4-1) |a'| 2 + |6'| 2 ' 



(51) 



In the entropy formula Eq.(45), we note that there is a degeneracy factor (d| — 1) associated with 
the \<j)' 2 ) state. Although in the example we have considered here, we took punctures carrying 
fundamentals, we see that the entanglement entropy is sensitive to the conformal block that the 
punctures in A (or B) may fuse to. This tells us if we label a state in terms of fusion or equivalently 
conformal blocks, we can read off the entanglement entropy directly from the representation around 
the interfaces. 



3.1.2 B with a* and a* 



To confirm this reasoning, let us compute carefully the second case, which consists of two funda- 
mental punctures in A. In this case, there are two states, shown in Fig. 14 The fusion rules 





i+i> 



i<|> 2 > 



Figure 14: The two states represented by Wilson lines connecting punctures. 

tell us that the representations cut by the interface are the symmetric and antisymmetric rank two 
tensor product of fundamental representations. Let's call them a and a), respectively. Of course in 
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P22 



Figure 15: Upon gluing to form pA, we find four matrix elements. Note that p\2 and p2i have the 
Wilson lines crossed in opposite senses. 



SU{2), they are the same as 9 and 0. As before, we first glue along B and get four possible density 
matrices, as shown in Fig. 15 In both pn and p22, we have two unlinked Wilson lines, so when 



they appear in trp^, they will not give rise to factors (as was the case for P22 in the last example). 
Now, pu and P21 have crossed Wilson lines with opposite orientation. In this case we have to be 
careful with the framing of the link and assign to each overcrossing a Dehn twist factor t = e 27Tlh& 
(see Appendix B). So, if a pi2 appears together with P21, they would just give an unlinked pair of 
lines, and thus the Dehn twists cancel each other. Thus only the difference in the number of each 
matters. To get the full tr pA n , we need to know the expectation value of a general braid with j 
crossings; we call this Xj. We already know Xo = So°(da) 2 and X\ — S^td^. If we denote by Li 
a pair of lines with i crossings, we have the skein relation aL+i + /3Lq + jL-i = 0, and then 



If we define q — e 

Xj + q^ + ^Xj-r = 

Using the notation [x] 
equation to get 



aXj + fiXj-! + 7^j-2 

27Ti/(N+k) we haye a 



0. 



(52) 



P 



1/2 dI1Q 



7 l/(2iV) 



7 l/2_ 



7 l/2_ 



-1/2 5 



thus 



_ i 
q 2" 



1 

" 2N 



2-^ — 2 — — (in which case d& 



7 l/2- 



-1/2 



5 n ° 



[N + 1][N] 



XT [2] 
This formula is valid for all integer j. Finally, 



+ (-<? 



1 — N 2 

[N]) and t q~^N 
^AN][N-l] 



' 2N I 



<? 

on 



^ jlkllUn 



n 



j-k-iy. 



(aa*) j (ab* 



[N][N + 1]. , i=k 
L JL -\a + bq2N 



2n + W]W 



' [2] 

<{a*b)\bb 
1] 



(d&n (53) 
we solve the difference 

(54) 



*\n-j 



-k-l 



Xi- 



q U 
<->0 



1+iV 



[2] ' * ' [2] 
We recognize in these expressions the quantum dimensions 

[N][N + 1] ,_[N][N 



a — bqiN 



2n 



da 



[2] 



[2] 



(55) 



(56) 
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The entropy then takes the form 

S A = 

where 

At = 



In So 



dco\\ InAi - daXz lnA 2 , 



(57) 



|a - &</ 2AT | Z 

do- 1 a + bq 2iv |^ + — bq 2iv |^ 

|a + bq 2N | z 
cfe|a + 2iv |^ + d^la - 2at |^ 



(58) 



Ai and A2 indicate an orthonormal basis, corresponding to the two conformal blocks. For the old 
basis, we have 



(<Pi\<Pj) = S °d a 
We can define a new basis as follows, 



da t 
t* d a 



W2) 



—q 2N Jd t 



fd&q~ 

/d^q^ q~2N^fd L 




(59) 



(60) 



Again we can calculate the fusion matrix in this case. ](/)[) and \4> f 2 ) are conformal blocks in 
the horizontal channel, while one choice of the conformal blocks in the vertical channel will be 

AT 



-|0i) and 



/S °d & lT±/ ~ y/S^d &y /d? & -l 

matrix for SU(N) k as 



(|0i) ~~ ^dalfa))- Using their relation, we can calculate the fusion 



a a 

a. a* 



dn 



^dn 



fdL: 



(61) 



When N — 2 there's no difference between a and a*, and it matches the result of the previous 
subsection. 



The wave function in the new conformal block basis is 




^x/d, 



In terms of wave functions under the conformal block states, 



Ai = 



1 



,/|2 



uf\2 



A 2 = 



d&\a'\ 2 + \V\ 2 ' d&\a'\ 2 + \V\ 2 ' 

giving the probability of finding a state in a given conformal block. 



(62) 



(63) 



3.2 S 2 with three quasiparticles 



There are many other cases that we could consider; generically, they cannot be represented by 
ordinary Wilson lines. The simplest such case is the three-punctured sphere; in terms of Wilson 
lines attached to the punctures, this would look like a 'string junction'. However, the density 
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matrix pa can be thought of in terms of Wilson lines. And of course given what we have learned 
in the previous section, we know that the entanglement entropy can be computed directly in the 
conformal block basis. 



For example, let's put d, a* and 9 on S 2 as in Fig. 16 Here the density matrix has the same 




Figure 16: Sphere with three punctures, chosen as representations a, a*, 6. 

conformal block as one of the states we found in the last section. Thus, up to a normalization 
factor, Pq a — —pi + d&p2- Following the same construction as in the last section, we will get 





Figure 17: The two density matrix elements in the case of the three-punctured sphere. 



^o = (-d & +d & r+(d & rd § , |= 

Clearly, there is a single conformal block contributing here. 



l-n 



*A(a,d*);6 



hi Sn 



(64) 



Similarly for A, A*, insertions, we find 



S A(a,a*);0 ~ m5 0°- 



(65) 



Given these examples and further thought about the general case, we can generalize the three 
quasiparticle sphere to 

S A(ij) ; k = m <S ° + In 4. (66) 



3.3 Does the topological entropy depend on the entire S-matrix? 



In the previous discussion we have found that the entanglement entropy depends on the quantum 
dimensions which are determined by the top row of the S-matrix. It is natural to ask if in other 
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computations the other matrix elements of S will also enter. We will argue here that the answer is 
essentially negative. As an example, let us consider a case in which one might have expected that 
the other matrix elements matter, the case of a torus in a non-trivial state i (e.g. a state created 
by a Wilson loop in representation Ri) with a puncture in representation Rj. One might anticipate 
that the entanglement entropy of the region represented in Fig. [9] would depend on . In fact it 
does not. The reason is that when computing trp^, any possible explicit dependence on Si 3 cancels 
out when properly normalizing pa- As a consequence the result depends on the 5-matrix only 
implicitly through the fusion numbers Nij k . Given the structure of this result, it appears that this 
is a general property. 

Another issue is the state dependence; in particular, we may have a situation in which the Hilbert 
spaces have dimensions greater than one, and hence there is at least implicit dependence on Nij k . 5 
A simple case of interest is four punctures on a sphere; suppose z, j are in A while fc, £ are in B and 
that both i x j and k x £ contain a block rn. If Nij m and N^ 171 are larger than one, 6 then pa is 
a matrix of rank min(A^j m , A^ m ). Each eigenvalue p a of pA contributes a factor — p a \n(p a /d m ) 
to the entanglement entropy (in the case of a single component interface). This would be summed 
over the possible fusion channels m. In the case where there are multiple interface components, 
we can organize the calculation into fusions A —> {mj} and B —> {mj}, where j = 1, . . . , / label 

the interface components. In this case, the rank of pA is min(A^^ m ^, A^ m ^), with N^™ 3 ^ the 
fusion number of punctures in the A region into the collection of blocks {rrij}. Each eigenvalue 
Pa of PA then contributes a factor — p a \n(p a / Ylj d mj ), which should be summed over a and the 
fusion channels. In this sense, the entanglement entropy depends on the fusion rules, but in an 
implicit way. This may be generalized to higher genus. In such a case, we also keep track of the 
representation along each handle, and these can make a contribution to the fusion numbers N; 
apart from taking this into account, the entanglement entropy is computed as we have described 
here. All of the examples considered explicitly in this paper (all of which had rank(pA) — 1) may 
also be expressed in this language. 



Chern-Simons Theory of the Topological Entanglement Entropy 
of Fractional Quantum Hall States 



The FQH states are topological fluids whose low energy effective field theory is a Chern-Simons 
gauge theory. As we saw in the preceding sections the entanglement properties of SU(N)k Chern- 
Simons gauge theories depend on the modular 5-matrix, which yield the quantum dimensions, and 
on the fusion rules of the excitations. We also saw that the entanglement entropy depends on the 
topology of the surface and on the regions that are being observed, and that when the states are 
degenerate the entanglement properties naturally also depends on which state is considered. 

In this section we apply the general results we derived in the preceding sections for Chern-Simons 
gauge theories to the computation of the entanglement entropies for both Abelian and non-Abelian 
fractional quantum Hall states. The results we will derive here apply only in the strict topological 
limit, that is for systems in the thermodynamic limit and for observed regions of size L much larger 



5 For SU(2) k , the Nij k 's are either zero or one, and so t his issue doesn't arise 



The explicit examples considered above in Section 3.1 have Nij™ < 1 
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than any intrinsic length scale of the physical system. In doing so we can only obtain the universal 
topological entropies. It will suffice to identify which Chern-Simons describes each case of interest 
and to use the results of the preceding section to compute the entropies. 



The entanglement entropy for the (Abelian) Laughlin FQH wave functions [41], as well as for the 
non-Abelian FQH pfaffian wave functions [121 EH], was calculated numerically recently in several 
papersjlU SS] which attempted to extract the topological entropy 7 for these states. This is in 
practice difficult to do numerically due to the large non-topological area term which needs to be 
subtracted. Similarly, the computation of the topological entropy in the conceptually much simpler 
Z2 topological phase of the quantum dimer model on a triangular lattice, which has a small but 
finite correlation length, presents similar difficulties [46J. The deconflned phases of 2+ 1-dimensional 
discrete gauge theories are actually the simplest models of topological phases [371 SHI EH]- For 
Kitaev's toric code state [23], i.e. the ultra-deconflned limit of a Z2 gauge theory, a state with a 
vanishing correlation length, it is simple to compute the entropy [50, 14J. The (non-topological) 
effects of a finite correlation length in a topological phase have been discussed in detail recently 
|51| . The scaling behavior of the entanglement entropy across a Z2 confinement-deconfinement 
phase transition was recently studied numerically [52J, as well as the role of thermal fluctuations 
on the behavior of the entropy in the Z2 topological state [53J. 

To proceed we will need to identify the Chern-Simons theory appropriate for the FQH state of 
interest. There is a well developed body of theory which does that and it is reviewed in Appendix 
A. The identifications that we need are the following: 

1. For the Abelian (Laughlin) FQH states, at filling factor v = 1/ra (with ra an odd integer) 
the effective field theory is an Abelian Chern-Simons gauge theory U(l) m (see Ref. jl2l 153]). 
It is straightforward to extend these results to the case of general Abelian FQH states. 

2. The bosonic non-Abelian FQH states are described by a Chern-Simons gauge theory for 
SU (2)fc, whereas the fermionic non-Abelian FQH states are described by Chern-Simons gauge 

theories whose CFTs are cosets of the form [SU@)/U(i)] 2 x U{1) (see Refe.02ll2ZlE3ll53l5Sl 
and Appendix A). 7 

3. The Chern-Simons theory describing generalizations of the p x + ip y superconductors (see 
Ref.dHlIinilSZ]) have a coset CFT (SU{2)/U(l)) k . 

4. The results presented here can be generalized to other non-Abelian FQH states of interest, 
e.g. the unpolarized non-Abelian states of Ref.[58] (and references therein) which involve 
more complicated systems such as £77(3)2 and others. We will not discuss these cases here. 



4.1 U(l) m Chern-Simons: The v = ^ FQH Laughlin states 

We will begin with a discussion of the v = 1/ra FQH Laughlin states which correspond to a £7(1) 
Chern-Simons theory at level ra. For the fermionic states ra is an odd integer, whereas for the 



7 This tensor product notation is ambiguous. The full RCFT has an extended chiral algebra. Its primaries are 
those in the tensor product which are local with respect to the current J + (defined below) with conformal dimension 



1 + M/2.[57l|43] In Section 4.3 we describe these structures in detail. 
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bosonic states m is an even integer. This case, and its connection with the modular 5-matrix and 
quantum dimensions, was discussed in great detail in Ref.[37j. For completeness, here we present 
only a summary of the relevant results. The description of the edge states of the Laughlin states 
in terms of a compactifled CFT is due to Wen. [59 , 60J. 

The U(l) theory consists of a compact free chiral boson of compactiflcation radius R. We will 
normalize 8 the field such that its correlator is ((/)(z)4>(0)) ~ —Inz. There is a U(l) current Jo ~ idcj) 
and operators Oq ~ exp(iQ0/ R) of conformal dimension Kq — Q 2 /2R 2 . If Q E Z, then Oq is 
single- valued. The characters of this model are 

q(n/R+wR/2) 2 /2 

Xn,w (r) = ^ , (67) 



where q = e 27TlT . In a rational CFT, we have that the radius is given by R — \j2p ! jp where p,p' 
are co-prime integers, in which case we can rewrite these characters as 

pp'{r+i/2pp') 2 

*,w = m . m 

In this expression, r E Z and —pp' < £ < pp' . These can be organized into characters of an extended 
algebra generated by Jo and operators J±, a set which closes under the action of the modular group. 
Generally, we find 

pp'{s+i/2pp') 2 

Xe(r) = J2 q V{T) ■ (69) 
where £ E —pp' + 1, ...,pp' .^T\ The modular 5-matrix for these characters is 

S/ = * (70 ) 

as can be easily established through Poisson resummation. We will use these formulae in later 
sections. 



In the case where pp' is even, this can be refined (that is, it is consistent (with respect to modular 
transformations) to consider a subsector of the Hilbert space) to 

q m(s+n/m) 2 /2 

X[n] (r) - Xi=2n + X£=2n+pp> = 2^ , (71) 

where we have identified (when pp' is even) m — pp' /2 and n is in the range 0, 1, m — 1. This is 
the case that obtains for the Abelian Laughlin states, 9 and we will refer to this theory as U(l) m . 
The set of primaries are in one-to-one correspondence with the states of a bulk Chern-Simons 
theory at level m. The extended current algebra is generated by J and J± ~ exp{±i^fm 0), the 
latter having dimension h± = m/2. J + is the operator that shifts n by m, leaving the character 
invariant and in the physical application, is interpreted as the electron. Requiring that primaries 
have local operator products with J±, we find 02n/p — exp(m0/y / m), of dimensions n 2 / (2m); these 
correspond to the fractionally charged quasiparticles. 

8 In string theory conventions, this corresponds to units a — 2. 

9 Specifically, we can take p = 2m, p = 1, which gives radius R — 1/^/m. These values are of course ambiguous 
up to T-duality, which acts as p <-> p , R — > 2/R. 
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Under a modular transformation, it is easy to establish (using Poisson resummation) that 

X[n<](-l/r) - E ^ e27rW/m X[n](r)- (72) 
Thus, we read off the modular 5-matrix 10 

and thus the total quantum dimensions 

v=(s °)- 1 

For the fermionic Laughlin states, m is an odd integer so that J± is a fermionic operator, which 
usually identified with the electron. In the special case m = 3 the electron has dimension 3/2 and 
the theory is equivalent to a superconformal field theory. For the bosonic Laughlin states m is 
even, and J± is a bosonic operator. In the special case m = 2, the theory is equivalent to SU(2) 1 . 




4.2 Coset [SU{2)/U(l)] k theories 

Chern-Simons theory whose CFT is the coset [SU(2)/U(1)]2 describes the two-dimensional time 
reversal breaking superconductors with symmetry p x +ip y . In some sense this is the simplest system 

with non-Abelian statistics. Here we will consider the general case of the coset [SU(2)/U(l)]k- 

To construct this coset, we begin with the SU(2) k characters, which are labeled by £ = 2j = 
0, 1, . . . , fc, and decompose them with respect to [SU(2)/U(l)]k x U(l) k 

k 

Xe U{2 \r)= J2 XeT^x^ir), (75) 

r=-k+l 

where Xr^\r) is given by Eq.([69]), with k — pp 1 ' . Since we know the modular transformations of 



both the SU(2) (Eq.Q) and 17(1) characters (Eq.([70])), we read off 

X^(-l/r)=5 ( , ir) ^) X S^(r) J (76) 

with 

I i r _/ n i -\\lnl i 1 M 



^^(1), r / 1 



k(k + 2) 



sm 



tt(£ +!)(£' + 1) 



fc + 2 



This should be restricted to r+£ = (mod 2) (since there is such a correlation between £7(1) charges 
and SU{2) representations). Note though that in this formula, equivalent characters appear twice, 



10n 



The <S-matrix of Eq.(82) is complex, symmetric and unitary. It differs from the result of Ref . [37] that found a 



real symmetric matrix. Nevertheless the quantum dimensions agree. 
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as xT-r et ~ X C k-lr±k- Accounting for these caveats, we find the independent characters, which lead 
to the final form of the 5-matrix: 



k(k + 2) 



sin 



7r(l+l)(l 7 + l) 



-iirrr' /k 



(78) 



For the case of most physical interest, we have k = 2, and the coset primaries may be taken to be 
(0;0), (1; 1) and (0;2). This is in fact just the chiral Ising model with (0; 0) ~ /, (1; 1) ~ a and 



r»coset 

^k=2 — 



IS 






1 


V2 


1 







-y/2 


1 


-y/2 


1 



(79) 



which agrees with the results of Ref . [37] . 



4.3 Moore-Read and Read-Rezayi FQH states: pfaffian and generalized parafermion 
states 



We now turn to the Moore-Read and Read-Rezayi non-Abelian FQH states, and their generaliza- 
tion. The filling factor of these states is v — k/(Mk + 2); M even corresponds to bosonic states 
and M odd to fermionic states [43J. As discussed above (and in Appendix A), these states are 
described by [SU{2)/U(l)} k x [7(1) CFTs, with a suitably defined level for the 17(1). Examples of 
these states are the well known Moore- Read pfaffian states. The the fermionic state with k = 2 and 
M = 1 has filling factor 1/2 (5/2 in the experiment), and the related bosonic state at filling factor 
v = 1 has k = 2 and M = 0. The states with k > 2 are the Read-Rezayi parafermionic states. 

We will discuss both the general fermionic and bosonic states with fixed k and M. The RCFT of 
interest is in all cases embedded in (SU(2)/U(l)) k x U(l) k ^ Mk+2 y We will consider the cases of 
k even and k odd as their structure is somewhat different. Here we only present details for the 
simpler cases. The details of the derivations for the general case are given in Appendix C. 

By reasoning similar to the above, the resulting 5-matrix can be obtained by multiplying coset 
and f7(l) characters. For the pfaffian state k = 2, the coset is a Z2 parafermion. The resulting 
5-matrix will, up to identifications, be given by 

S^sZ'^ = (Sr et )(e,r) {i ' y \S U{1) * M+ *)/. (80) 

Primaries of this theory will be given by products of the Z2 primaries {/, <r, ip} with ?7(1)4m+4 
primaries of the form 0^/ p . We seek a set of such operators that close under operator products 
and are loca l wit h respect to a suitable extended current algebra, which will be generated by 
J± rsj ijj e ±WM+icf>^ w here ^ is the Maj orana fermion of Z2. For simplicity, we will consider two 
cases here, M = (take p = 2, p f = 1, radius R = 1) and M = 1 (take p = 4, p' — 1, radius 
R=^lj2). 

In the case of k — 2 and M — 0, we find the integer-weight J± ~ ^ e ±2( ^ as suitable extended 
currents. Requiring locality of operator products, we then find that the primaries of this theory 
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are given by /, ^, cre^/ 2 (all others are related to these by action of J±). These in fact are just 

the primaries of S77(2) 2 , as we should expect. This is the bosonic pfafflan state. The associated 
modular 5-matrix was given in Section II and Appendix C. 

In the case of k = 2 and M — 1 (the fermionic pfafflan state), we find J± ~ i\) e ±l ^^ as suitable 
extended currents. Requiring locality of operator products, we then find that the primaries of this 
theory are given by 

/^tre^/^e^/A (81) 

This set closes under fusion (up to the action of J±). These operators have weights 11 (0,0; 0), 
(0, 2; 0), (1,1; ±1) and (0, 0; ±2) respectively. We can then read off the <S-matrix: 



2y/2 



( 1 
1 

V2 
V2 
1 

\ 1 



1 
1 

-V2 

1 
1 



V2 
-V2 





-iy/2 



y/2 
-y/2 



+iV2 



1 
1 

+iy/2 
-iy/2 

-1 

-1 



1 
1 

-iy/2 
+iy/2 

-1 

-1 



(82) 



from which one can read-off the total quantum dimension is V — 2\/2. This model can also be 
viewed as (the NS sector of) an N — 2 superconformal current algebra, as the current J+, having 



conformal weight 3/2, can be viewed as a supercharge. So the operators listed in Eq.(81) are then 
viewed as superconformal primaries. [^3} 157] 



We will now consider the interesting example of the parafermionic states at k — 3 and M = 1: 
the Read-Rezayi parafermionic state for fermions at filling factor 2 + 2/5. The k = 3 coset has 
primaries at (£, r) = (0, 0), (1, ±1), (2, 0), (3, ±1), which we will refer to as /, <j±, e, rp± respectively. 
Explicitly, denoting s p = sin(7rp/5), we have 





( 81 


S2 


S2 


S2 


si 




si 






S2 


e -i7T/3 5i 


e +i7T/3 5i 


-si 


_ e -i7r/3 S2 


—e 




2 


S2 


e +i7T/3 5i 


e -i7T/3 Si 


-si 


_ e +i7T/3 S2 


—e 


-i7T/3 S2 




VIE 


S2 


-si 


-si 


-si 


S2 




S2 






Si 


_ e -i7r/3 S2 


_ e +i7T/3 S2 


S2 


_ e -i7r/3 Si 


—e 


fi7r/3 s 






\ si 


_ e +i7T/S S2 


- e -™/3 S2 


S2 


_ e +i7T/3 Si 


—e 


-m/3 Sl 


) 



(83) 



For this case there is an extended algebra generated by the h = 3/2 operator Q+ = ^+ e ^<P/Vi5^ 
where is a free boson of the U(l) theory that we are attaching |57j . Representative primaries 
are (e,r;s) = (0,0; 0), (3,-l;l), (3,1; 2), (0,0; 3), (3,-1; 4) and (*,r;a) = (2,0;0), (1,-1; 1), 
(1,1; 2), (2,0; 3), (1,-1; 4). One can check that these have local OPE's with Q+ and are closed 
under fusion. As we will see, it is convenient to group them into groups of k + 2 = 5, as 
given. The theory obtained this way is actually an N = 2 superconformal theory, with super- 
charges Q± (Q_ being ^_ e ~^l^). Q + groups collections of conformal primaries together, i.e., 
{(0, 0; 0), (3, 1; 5), (3, -1; 10)}, {(3, -1; 1), (0, 0; 6), (3, 1; 11)}, {(3, 1; 2), (3, -1; 7), (0, 0; 12)}, 
{(0, 0; 3), (3, 1; 8), (3, -1; 13)}, and {(3, -1; 4), (0, 0; 9), (3, 1; 14)} and {(2, 0; 0), (1, 1; 5), (1, -1; 10)}, 
{(1,-1;1),(2,0;6),(1,1;11)}, {(1, 1; 2), (1, -1; 7), (2, 0; 12)}, {(2, 0; 3), (1, 1; 8), (1, -1; 13)}, and 
{(1, —1; 4), (0, 0; 9), (1, 1; 14)}. Each of these triplets represents a superconformal family. When we 

n The notation (£,r;s) represent the coset weights (£, r) and the £/(l)-charge s. 
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compute the 5-matrix with respect to the extended symmetry, we treat these groupings as one. 
That is, computing the 5-matrix element on the grouping gives a 3 x 3 identity matrix times a 
factor. We collect those factors into the following 5-matrix. 



nFRR 
C> k=S 



2 / sin(7r/5) sin(27r/5) 
5 \ sin(27r/5) — sin(7r/5) 



/ii i i i \ 

1 CJ2 ^4 ^1 ^3 

1 UJ4 UJs (j02 W\ 

1 UJ\ U2 (j0% UO4 

\ 1 CJs UJ\ CJ4 CJ2 ) 



(84) 



where we have used the U(l) 5-matrix is S s s ' = —L^e 27 ™ 55 '/ 15 . Above we used the notation is 



u p = e 27Tip / 5 . The coefficient out front is ^= • ^= • the factors being the coefficients of the coset 
5-matrix, the U(l) 5-matrix and the order of the automorphism (5 in 15), respectively. Note that 
it is easy to read off then the total quantum dimension 



V 



= V=^5 + 5(s 2 /si) 2 = v/5(i + 2 ) 



5 ° 2sin(7r/5) 

where here (p = (y/5 + l)/2 denotes the Golden Ratio (not the chiral boson!). 



(85) 



In Appendix C it is shown that for general k and M, the primaries are the highest weight states of 
the form 



%/-2[ft]) exp 



.l + nfc- {Mk + 2)[fj} 
^k{Mk + 2) 



,£ + nk- (Mfc + 2)[^ 11 



<t> or ^_ 2[ ^i D exp i 



y/k(Mk + 2) 



(86) 

where ip^ r ) are Z^-parafermion primaries, n and £ are integers (with a suitable range, see Appendix 
C), [x] is the closest integer to x. For general k the 5-matrix is given by 



^(k + 2)(Mk + 2) 



sin 



^(l+l)(f + l) 
fc + 2 



exp 



m(-MW + 2ln ! + 2f n + 2W) 



Mk + 2 



One can read off from this the total quantum dimension V for all M and fc, since 



(87) 



1 

V 



<^>0 



sm 



^(fc + 2)(Mfc + 2) ' \k + 2 



7T 



In Appendix C we also show that for general (odd) k, the modular 5-matrix may be put into the 
simpler form 



o {l';n'} 



V(A: + 2)(MA: + 2) 



sm 



^(l+l)(f + l) 
k + 2 



(-if exp 



2niknn' 
Mk+2) ' 



(89) 



with ^ = 0, . . . , and n = 0, . . . , M/c + 1. The result of the fc = 3, M = 1 example considered 



above matches with this. 
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5 Conclusions 



In this paper we computed the entanglement entropy for Chern- Simons gauge theory with general 
gauge group. We considered the specific cases of SU(N)k and various cosets (of interest in the 
theory of FQH states). We have done this by direct computation in the 2 + 1-dimensional gauge 
theory using surgery techniques. We found that the entanglement entropy for these theories can be 
written as a Chern-Simons path integral in a complicated 3-manifold whose details depend upon 
the topology of the spatial surface and the way it is partitioned to compute the entanglement. In 
all cases the entanglement entropy can be expressed in terms of modular features of the dual two- 
dimensional conformal field theory. We found that in general the entanglement entropy depends on 
the universal data of this topological field theory (that is the quantum dimensions, the fusion rules 
and the corresponding fusion numbers). However in cases in which, due to topology the ground 
state is degenerate, the entanglement entropy may also depend on the choice of state. 



Acknowledgments 

We thank Paul Fendley, Eun-Ah Kim, Alexei Kitaev, Michael Levin, Joel Moore, Chetan Nayak, 
Stefanos Papanikolaou, John Preskill, and Kirill Shtengel for many discussions. This work was 
supported in part by the National Science Foundation through the grant NSF DMR 0442537 (EF) 
and the US Department of Energy grant DE-FG02-91ER40709, and the Stanford Institute for 
Theoretical Physics (EF). 



Appendix A: Chern-Simons gauge theory and Fractional Quantum 
Hall States 



In this appendix we review the connection between FQH syes and Chern-Simons thory. The FQH 
states are incompressible electron fluids which, due to the presence of the large magnetic field, 
have an explicitly broken time reversal invariance. These topological fluids have a ground state 
degeneracy which depends only on the topology of the surface on which the fluids reside. [62] In 
their low energy and long distance (hydrodynamic) regime, the FQH fluids behave as topological 
fluids. The excitations of these fluid states ( "quasiparticles" ) are vortices which, in general carry 
fractional charge and fractional (braid) statistics. [63] In other words, the long-distance correlations 
in this topological fluid are encoded in the fractional charge of its vortices and, more significantly, in 
the non-local effects of fractional statistics. However, since these topological fluids are condensates 
of electrons, these quasiparticle (vortex) states are local with respect to states representing an 
electron. This condition, and the quantization of the electron charge e play a key role in the 
properties of the effective theories of these fluids [42J. This structure is also responsible for the 
extended symmetries in the corresponding CFTs that we described in the body of the paper. 

In the hydrodynamic regime i.e., at energies low compared to the quasiparticle excitation energies 
and on length scales long compared with the magnetic length, the physical properties of these FQH 
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fluids have been shown [6U [65l |66j [671 EH EH1 E3 EE] to be described by an effective (topological) 
quantum field theory, the Chern-Simons gauge theory in 2 + 1 dimensions. [9| 169] 



A.l The Abelian Quantum Hall States 

We will consider first the Laughlin states, whose wave functions for a system of N particles at 
filling factor v = 1/ra are [H] 

N 

-£m 2 /^ 2 

V m (zi,...,z N ) = '[[(zi-z j ) m e <=i (A-l) 

i<j 

where {zi} are the complex coordinates of N particles, and t is the magnetic length. Following 
Ref.[54], we write the effective field theory of the Laughlin FQH states, which have filling fraction 
v = 1/ra (with ra and odd integer for fermions and an even integer for bosons), as U(l) m Chern- 
Simons theory, whose Lagrangian density is 

772 

C = —e^xA^Ax (A-2) 

up to irrelevant operators whose effects are negligible in this extreme infrared regime. The field 
strength of the gauge field A^ is essentially the hydrodynamic charge current of the topological 
fluid 

J» = -^v\dvA x (A-3) 

where e is the electric charge. The (infinitely) massive quasiparticle bulk excitations, the vortices 
of the topological fluid, are represented in this limit by temporal Wilson loops (the world lines of 
these quasiparticles) . On a manifold of genus g, the ground states of the FQH fluids are degenerate 
[62]. In the Chern-Simons description, the degeneracy is ra 5 ', where ra is the (quantized) level of 
the Chern-Simons theory. 

This hydrodynamic description generalizes to describe all other Abelian FQH states |70l[71]. The 
corresponding effective field theory is a Chern-Simons gauge theory of a tensor product of U(l) 
gauge groups at various levels, as well as non- Abelian groups at level 1 (which only have Abelian 
representations of the braid group) [68, 72]. Here we will consider only the simpler Laughlin states 
as the generalizations of our results to the other Abelian fluids is straightforward. 

A direct consequence of the topological nature of the FQH fluids is that, for a physical system with 
a boundary, their excitations are gapless. These edge states are described by a chiral conformal 
field theory (CFT). [S3 [59l EH] In particular, there exists a one-to-one correspondence between the 
gapped bulk quasiparticles and the primary fields of the edge chiral CFT. At the edge, the effective 
CFT for a Laughlin state is a U(l) m chiral boson (j) in 1 + 1 dimensions. 

The physical requirement that the quasiparticle edge states are local with respect to the electron 
operator leads to the compactiflcation of this CFT, which now becomes a chiral rational conformal 
field theory (RCFT).[42l [54] Thus, the only operators allowed in the edge chiral RCFT must obey 
the condition of being invariant under the compactiflcation condition 0^0 + 2nR, where R is the 
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compactification radius (see Section [I]). The compactification condition leads to a truncation of 
the spectrum which has m distinct sectors, the same as the ground state degeneracy on the torus. 
Since in the Abelian states all the states are one-dimensional representations of the braid group, 
it follows that all the quantum dimensions of these states are equal to unity, di = 1, where here 
i = 1, . . . , m. 



A. 2 The non-Abelian Quantum Hall States 



The non-Abelian FQH states are more interesting and have a more intricate structure. Although 
no simple unified effective field theory of all the non-Abelian states yet exists, in all the cases that 
have so far been studied the effective field theory contains at least a U(l) m gauge group at some 
level m (for the charge sector), and a non-Abelian gauge group such as SU(2) k at some level k. 
For instance, for the Moore-Read pfaffian FQH states [32] 

N 

* q (z u • • • , z N ) = Pf l—^-) II (* - ^ e 1=1 ( A - 4 ) 

x J ' i<j 

where Pf (^ Z .^ Z ^J is the pfaffian of the matrix. It has long been known [42J that this wave function 
can be regarded as a correlator in an Euclidean two-dimensional CFT. Indeed, the Laughlin factor 
is simply the expectation value of a product of vertex operators of a chiral Euclidean boson = 
e iy[rn$ ^ m a neutralizing background) and a correlator of majorana fermions ij; in an Ising chiral 
Euclidean CFT. The theory of the edge states of the Moore- Read states is a chiral CFT described 
in Section [U 



It turns out that for the bosonic state at m = 1 (with filling factor v = 1, which may hopefully be 
accessible to experiments in ultra-cold gases of bosons in rotating traps [20J), the effective field theory 
in the bulk is simply an SU{2) 2 Chern-Simons gauge theory [271155], without any U(l) factors. The 
excitation spectrum of this state consists of a particle with non-Abelian braid statistics, the Moore- 
Read "non-abelion" , created by the Ising primary field cr, whose quantum dimension is d a = \/2 
(since the degeneracy of a state with 2n such vortices is 2 n ~ 1 [73]), a Majorana fermion (with 
quantum dimension = 1), and the identity field / (the boson also with quantum dimension 
do = l)- 

The two-dimensional time-reversal breaking superconducting state with symmetry p x + ip y , appar- 
ently observed in Sr2Ru04, is also closely related to paired-Hall pfaffian states [75J. It is also a 
non-Abelian state and it is also, up to an U(l) factor, an SU{2) 2 state. Indeed, Fendley, Fisher and 
Nayak [37] have computed the effective quantum dimension for this case as well and found that it 
is also equal to \/3. 



The fermionic pfaffian state at m = 2 is the natural state to explain the observed plateau in the 
quantum Hall conductance at filling factor ^ = 5/2 = 2 + 1/2. As can be seen from the structure 
of the Moore- Read states, c.f. Eq.(A-4), the bosonic state at m = 1 and the fermionic state 
at m = 2 differ only by the Laughlin factor. Hence, one expects the fermionic state also to be 
connected to SU(2) 2 which, up to some caveats [27], is essentially correct. The SU(2) symmetry 
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of the bosonic case is dynamical and it is broken in the fermionic case due to the change in the 
compactification radius of the boson (see Refs. [271 E23 EE]-) Thus, one expects the the quantum 
dimensions of SU{2) 2 should play a role here too. However, the presence of the additional U{1) 
factors, associated with the charge sector, and the breaking of the SU (2) symmetry changes the 
dimensions. Recently, Fendley, Fisher and Nayak [37] have analyzed this case in detail. The upshot 
of their analysis is that the fermionic pfafflan state has a total of six primary fields: the identity 
/, two (conjugate) non-abelion primaries ae ±lc ^^ 2y ^\ the Majorana fermion -0, and the Laughlin 
quasiparticle and quasihole e ±l ^^^. The quantum dimensions of these states are, respectively, 
dj = 1, d(j — v^2, = 1, and 1 for the Laughlin vortices. 

More interesting from the point of view of topological quantum computing, but not yet clearly seen 
in quantum Hall experiments, are the Read-Rezayi parafermionic states. [13] The Read-rezayi states 
are constructed in a manner analogous to that of the Moore-Read states. The main and important 
difference is that the pfafflan factor, which as we saw is equivalent to a correlator of Majorana 
fermions in a chiral Ising CFT, is replaced by a parafermion correlator of a parafermionic chiral 
CFT. In particular, the simplest bosonic parafermionic state can be represented in terms of the 
Chern-Simons theory SU(2) 3 [55]. The fermionic counterpart can also be understood in similar 
ways. The interest in this state stems from its non-Abelian vortex. As a consequence of their fusion 
rules, the topological degeneracies of these vortex states follow the Fibonacci sequence, and the 
quantum dimension of this non-Abelian vortex is the Golden Mean [25] . 



Appendix B: Calculation of the SU (N)k modular <S-matrix and Fram- 
ing Factor 



B.l The 5-matrix Elements 



We use (14.247) in [EE] 



^6 



Xx 



-2ttz(5 + p) 



(B-l) 



to compute the 5-matrix elements. The notations used in this Appendix are also adopted from [6T] 
and differ somewhat from the notations we used in the body of the paper. 



For SU(N) k , there is a natural orthonormal basis for the root lattice to compute the characters. 
It is constructed as follows. Pick N dimensional unit lattice with unit vectors {e^}, i = 1, . . . , TV, 
then the simple roots of SU(N) can be written as 



Oii €i 



i = l,. 



(B-2) 



i.e. the root space sits in the N — 1 dimensional subspace with J2iLi n i ~ f° r an y element 
^2iLi n i e i- Any integral representation A = Ya^i* ^ expressed in the orthonormal basis, 

becomes 



N 



(B-3) 
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where k = X^jL/^j * s the partition of the associate Young tableau, and k — ~^Yl!j=i 3^j- I n 
particular, the Weyl vector becomes 



N-l N 



p=T, u * = - i - ^te- ( B - 4 ) 



Since 



Xx 



-2ttz(5 + p) 



D 



A+p 



k+9 J _ J2 w ew e M e 



w(\+p),( 



-2ivi(5+p) 
k+g 



-27Ti(5+p) \ 

k+g ) 



[wp,( 



k+g 



(B-5) 



— 2-Ki 

where Weyl group W is simply the symmetric group of {e^}, if we define q = e k +a , we will have 



^5 



E 



N „(AT_ s ._^ p )(^ (5 . + Ar-i-^ (5 -^ p ) 



e^ e ( 5 ) rui ? 

det[g^ A * +Ar_i_ ^ A_ ^ p ^ /(5 ^ +Ar_J ' _ ^ <5_ ^ p ^] 



= q 



c ] Le ^q( N - i -Kp)(l5j+N-j-K5-Kp)j 
-« A Sf =l(t ,, + *-j-,«-^ 

det[q( N ~ i ^ ls i +N ~ J ~ K5 ~ Kp>) ] 



det[q( lxi+N ~ i ^ l5 i +N ~i~ K5 ~ Kp " ) ] 

det[q^ N ~ i ^ ls 3 +N ~^~ Ks ~ K/ p^] 
q f-^-^^ lxj det[ g (^+^fai+^-i)] 
9 det^"^^^')] 
q-N(K S +K p )K X g x ({q l5 i+ N ~i}), 



where in the last step we have used the Schur function 



det \x n 



det \x 



N-ii 



(B-6) 



To calculate the Schur function specialized at {q ls i +N J }, we can use one of the Giambelli's Formula 
(see Appendix A.l around (A. 5) of [77J for detailed discussion), 



S x ({xj}) = det[£? A r, 



(B-7) 



where A T is the transposed partition of A, and are the elementary symmetric polynomials 
generated by 



N 



E(t) = H(l + x i t) = ^2E m t m , 



(B-S 



2=1 



m=0 



and £y = for j < 0. 
Now, in our case 



N 



E(t) = H(l + q 



(B-9) 
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and in principle by expanding E(t) we can read off all the E m and calculate S\({q ls i +N J }). From 
(14.217) of [EE] we can see that <%° = thus 



^_6o_ 
c o c o 



Sc — St 



(B-10) 



To calculate S °, we need to use (14.217) of [61J. For SU(N) k , |A+| = M^zii, |p/Q-| 
r = iV — 1 , thus 



N,g = N, 



jy(jy-i) 

= i 2 



. / a - NiK " )2 det\a iN - i)iN - j) ] 

\jN(k + N)^ q aet[q 1 



N(N-l) 

i 2 



N(k + N) 



N-l l 



1 



l<i<j<N 



l<i<j<N 
sm jt)- 



The quantum dimensions are as follows, 



rf A = ?0 = q- NK ^S x ({q N -i}). 



First we calculate the generating function, 

N N m n r _i 

Eit) = n(i + q N -H) = i + y: t m n 

z=l m=l r=l 



(B-ll) 



(B-12) 
(B-13) 



and then plug into 

Sx({q N - j })=det[Eyr i+j _ i \. 

For example, for the fundamental representation a of SU(N), the transpose a T = {1, 0, 0, ... , 0} 
and 

/Ei E 2 En \ 

E Ei ■■■ 
E ■■■ =Ei 



S Q ({q N -i}) = det 



\ 



q N -l 



q-1 



(B-14) 



E J 



Thus 



- N *-^q N -l q N / 2 -q~ N / 2 



(B-15) 



cc/2 —x/2 

Here we have used the g- number notation, defined as [x] = ^ 1/2 ~^_ 1/2 . Using this, we can write 



as 



m q N - q 1 - 1 _ {\ n=i [N + 1 - r] _ (N-i)m ™ [jV + 1 - r] 
g'-l " : " [r] _<? 2 11 M 

r=l r=l L J r=l 1 J 



^=nv^r- = n< 



(B-16) 
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We also want to check d^*, for which a* T — {N — 1, 0, 0, . . . , 0, 0}, and 

/ E N _ X E N • • • \ 
Eq Ei 
S a *(U N ~ j }) = det E ••• 

£ / 



V 



(Af-l) 2 

= ^iV_l =</— a— [JV]. 



(B-17) 



Thus 



d a . = — ^+^^[AT] = [AT] = d & . 



(B-18) 



For symmetric and antisymmetric rank two representations a and Co we have a T = {1, 1, 0, . . . , 0, 0} 
and uj t = {2, 0, 0, ... , 0, 0}, thus 



S a ({q N - j }) = det 



( E\ E2 
Eq Ei E2 
E 



En \ 



v 



E(-E 2 = q N -\[N] 



[2] 



and 



Mti N - j }) = det 



••• E J 

( E2 E3 ■ ■ ■ 

E E x 

E 



En \ 



V 



= E 2 = q 



[2] 



E J 



Thus 



and 



da — q 2 N 



+ _ [N][N + 1] 



[2] 



[2] 



, _ N n=i i +(N ^ [N][N - 1] _ [TV] [TV - 1] 

[2] " [2] ' 



[2] 

(B-19) 
(B-20) 



(B-21) 



(B-22) 



In the weak coupling limit when k — ► oo, [x] — * x, the quantum dimensions we just calculated 
match their classical values. 



The last thing we want to calculate is Sa a - Following the procedure described before, the partitions 
a = a T = {1, 0, 0, . . . , 0}, so 



S a ({q l «i +N -j})=E 1 ({q l «i+ N -j}), 
where £i({^ +jv ^}) = J2?=i q la ' +N ~ j = q N + = q N + g^[iV - 1]. Thus 



(B-23) 



C ol Q ot Q 
<? C o " fi 



o - c = Q- [1+ ^^(q N + q^W -imN} = q -^-H q ^ + [N -1})[N}. (B-24) 
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B.2 The Framing Factor and the Skein Relation 



For any representation A, the unit of Dehn twist factor is t = e 2nih x, where 

(A, A + 2p) 



^ " %T#' (B - 25) 

For the fundamental representation a of SU(N) k , 

* 2(fc + c?) 2(fc + iV) 2iV(fc + A^)' 1 > 

where F here is the quadratic form matrix of SU{N) 1 which is the inverse of the Cartan matrix, 
and 

1-iV 2 

t = q-2ir. (B-27) 

In order to solve aL+i + (3Lq + 7L-1 = 0, we complete the path integral in two different ways and 
get 

^td & + (d a ) 2 + lt*d & = 0, (B-28) 

?(4&) 2 + **^ + ^ = °- (B-29) 
P Poo 



and 



Notice unlike in [9J, we keep the framing factor explicitly here. 

« = f 1 -**? = f l-W + g-^[N- 1] [N + 1] [N] (g 1 / 2 - q~^) = g~± 

P t 2 ^-d & -q-T[N-i\[N + i\(qV*-q-V2) gi/2_ g -i/2' 



7 = (d a ) 2 - 1 = [iV] 2 - 1 = = ( a m , 



Appendix C: 5-matrix of the (sU(2)/U(lj) x t^T) fc(Mfc+2) RCFT 

In this appendix we will calculate the primaries and 5-matrix of the (^SU(2) /U(l)^j x C/ (l)fc(Affe+2) 
theory. For simplicity we will assume k is odd for the moment, and discuss even k later on. The 
fields in the Z^-parafermion CFT, (the coset (SU(2)/U(1)^J ), are labeled by the SU{2) charge 

and its U(l) subgroup charge (£, r), where we take them as twice the traditional values, so that 
they are integers. For SU{2)^ 1 £ = 0, . . . , k and I — r = mod 2. We also have the identification 

(£, r) = (k - £, r ± k) = (£, r ± 2k) . . . . (C-l) 

Using this identification we can always map (£, r) for i > to (k — £,r ±k). So we can restrict to 
< £ < *±±. We will use ip(t^) for the corresponding fields and X(£,r) their characters in the coset 
theory. 
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With these properties, one can see that the states form Z^-loops generated by ij) or ^t, where ij) 
and ^ are the parafermion fields that appeared in the original SU(2) k currents J + ~ ^^VV^ 
and J_ - ^t e -i«v^A. We can identify -0 ~ 0(0,2) an d ift ~ ^(0,-2) f° r which /i^ = /i^t = 1 — 
The other fields in the SU(2) k , ^^ r) e ir ^^, have [7(1) charge r. The whole multiplet under the 
SU(2) k current algebra has character 

-.00 -j k—l 00 

E xw-V- 2 "*' 4 * = ^ E E <,"- 2 ' + « 2 /« (c-2) 

/V ' £=-oc /V ' t=0 p=-oo 

We have used the periodicity of X(£,r)- Also, the fields tfj(£^e l ( r+k ^/^^ satisfy the same locality 

condition as ^^e^/^^, thus they are in the SU(2) k theory as well. In fact, they are the k — £ 
multiplet with character 

-j k — l 00 -j k—l 00 

iEw E /- 2m+2p/£)2/4fc -^Ex^-^-^) E /- fc+2t+2 ^ 2 / 4fc . (c-3) 

^ ' t=0 p=-oo ^ ' t=0 p=-oo 



The highest weight state in V^r) with r = — £, . . . , £ and ^ = 0, . . . , k has h — 4^+2) ~~ wn i cn 
means for any integer p the highest weight state in ip(e i e+2 P k-2t) nas weight 

( £(£+2) (£-2t) 2 .r ._ n ; 

r _l 4(^+2) ~ 4fc 11 r-U,...,A 

^,/+2^-2t) - < (fc-l)( fc -l +2 ) _ (g-gt+fc^ : f , _ / , 1 u _ i ^~ 4j 

I 4(fc+2) 4/c 11 ^ — « "I" 1 



From the modular transformation property of SU(2) k and J7(l) we can also get that of the 
parafermions as follows 

W' V) = , 2 sin ^ + / )( ^ + 1) e-^/ fc . (C-5) 

Now, instead of J±, let's use Jj_ M/2 ~ ^eV 2 / fe + M * and ~ ^t e -V 2 / fe + M * The subscript 

indicates the weight of the currents. The same locality condition tells us that V>(^ r ) has to be 



multiplied by exp 



--(kn + r) 



for any integer n. Together with the Jo ~ id(/) and other 



fields with £ = 0, the new currents form an extended chiral algebra. The character under this 
symmetry is now 

I ^Z^ [l-(Mfc+2)t+fcn+pfc(Mfc + 2)] 2 

E 1 2k(Mk+2) ■ (°- 6 ) 

^ ' t=0 p=-oo 

The independent multiplets correspond to n = 0, . . . , M/c + 1. When M = 0, we get the original 
SU(2) k theory, as above. When M ^ 0, a little calculation shows that if < n < M£, the highest 

, l+nfc-(Mfc+2)[^] ^ 



weight state in the multiplet is V^-2[^])e ^ fc(MH2) , while if M£ + 1 < n < Mk + 1, the 



^ +nfc _(Mfc + 2)[^] 

highest weight state in the multiplet is ^n-i^e v /fc ( Mfc + 2 ) ? where [x] denotes the closest 
integer to x. We will call this character X[£,r(£,n);s(£,n)]i where 

r(£n) - { if 0<n<Ml 

l ' ) ~\i-2[^} if M£ + l<n<Mk + l 
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and 

sUn) = i t + + if 0<n<M£ 

{ > ] \ £ + nk- (Mfc + 2)[^] if M£+l<n<Mk + l [ ' 

are the corresponding U(l) charges of the primary fields. 

Under modular S transformation, the standard Poisson resummation result tells us that 

1 oo 9 N— 1 oo 9 

1 (r+piV) v-^ 1 27rirs 1 (s+piV)^ 

-— - > g 2iv _> \ __ e at — — \ q 2N (C-9) 

/v J p=-oo s=0 v /v J p=-oo 

which gives the S transformation property of the combined theory as follows, 

X[£,r(£,n);s(£,n)](-^/ T ) 

' ^ A ' ' 2 . n (£ + !)(£' + 1) 

-. sm 



v y v 

fe^o^o V(fe + 2)(Mfc + 2) ^ + 2 

fc(M/c+2)-l , , / x oo 

27ri(£-2t-Mkt + kn)s Txi(t-2t)(t' -2t' ) \ ( pi ff_Off\[T) 
\ P fc(Mfc + 2) fc H ' ; V 7 \ ^ 2fc(Mfc + 2) 

Zv ?? ( r ) Z. 9 

,s=0 v y p=— oo 

fc-l-i 

^M/c+l 7r( ^ +1)(f + 1) 
sm 



EE 



v '(* + 2)(Aft + 2) f f „'e * + 2 

27ri(l+fcn)(l / -2t / +n / fc) nU^' -2t') y^; />/_ 0f n 

\ p fc(Mfc + 2) fc /v^ ±i j x 2fc(Mfc+2) 

Z^ e ^( r ) 9 

t'=0 ' v ' p=-oo 

k+l 

2Mk+l 7r ^ +1 \/^/ + 1 \ T , ( -Mtt' + 2fa'+2<'n+2 t nn') 

2^ 8111 6 Mfc+2 X[f, r (f,,'); S (^,n')]( T ) 



#+1^+2)^0,^0 fc + 2 



i.e. 



S^.^ ' n j = , sm — -e Mk+2 . (C-10) 

V(fc + 2)(Mfc + 2) A: + 2 v 



We have repeated used the periodicity of x(^ r )- To make the formula simpler, let's define 
show that 



1 ^ 2 1; , thus since A: is odd, £ + cr^fc is always even. Define n — > n — — ^+ a ^ M ; then we can 



5(^ ;n )^ ,n j = / /7 sin — -(-1)^ e . (C-ll) 



7r(l+l)(l' + l), 
^(fc + 2)(Mfc + 2) "~ fc + 2 
From the periodicity, we can still have n = 0, . . . , Mfc + 1. Notice if we use the new parameter, the 
function of U(l) charges of highest weight states in terms of n has to be modified. The 5-matrix 
is now factorized into £ and n parts. For fc = 3 and M — 1 for example, this result coincides with 
that given in the text, eq. (84). 

When fc is even, the loop with £ — fc/2 has only fc/2 elements. The character of the multiplet in 

(kn+k/2)<t> 



which the state ip(k/2,k/2) e v /fe ( Mfc + 2 ) lives is 



k 

I 2 ' ^ [g-(Mfc+2)t+fcn+p "^^^' ; ] z j_ At_J - ^ [§ -(Mfc + 2)t+fcn+pfc(Mfc+2)]^ 



2 ~ L OO rfc / a /f i ou i ^ i ^ fc(Mfc + 2) l2 k-1 



„ (T) E E 9 2fc(Mfc+2) =^yEx ( |,|- 2t ) E ^ 2fc(Mfc+2) 



(C-12) 
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for n = 0, . . . , Mk/2. Apart from the range diffence, the character looks just the same as before. 
So the calculation of the modular 5-matrix is almost the same as the odd k case. 



X[£;r(£,n);s(£,n)}(-1/ T ) 

^ ^ S(£',n)^ ^ ) X[^,r(^,n / );5(^ / ,n / )]( T ) + 
Z'=0 ra'=0 



~2~ M/c+1 



v /(fc + 2)(Mfc + 2) v ? 5 



(2^+ ) sm " 



=0 „'=M*l + i 



fc + 2 



- — 1 / \ 

7T2(-M£|+2£n / + fen+2fcnn / ) 2 fe_2t / )V T ) [| -2t / -Mfct / +n / fc+pfc(Mfc+2)] 2 



Mfc + 2 



£ 

t'=0 



»7(r) 



2fc(Mfc+2) 



P= — (X) 



^ ^ S{i-n)^ ^ ) X[^,r(^ / ,n / );s(^ / ,n / )]( T ) + 
Z'=0 n'=0 



Mfc 

2 i , / 1^ 



V(fc + 2)(Mfc + 2) n ^ 2 



iri(-Ml^+2ln' + krL+2knn') 
> Mfc + 2 



fc-1 



t'=0 



X(k fc_2t / )( T ) °° [|-2t / -Mfct / +n / fc+pfc(Mfc + 2)] 2 



7?(t) 



E 

p=— OO 



2fc(Mfc+2) 



Mfc 
2 



S 5 (^;n) ( ^ ;n,) X^,r(^,n'); S («',n')]( r )+ ^ S^^'^^^k^^k^r) 



(C-13) 



^=0 n'=0 



n'=0 



We saw that the 5-matrix is the same as k odd case. The only difference is that when k is even, 
the £ blocks are not of the same size. To simplify, the best we can do is to make the redefinition 
n^n- £M/2, then 



c WW) _ 



. ir(£ +!)(£' + !) 



sin 



i M££' ( _ yni'+in' g nv7fc+^ . (C- 14) 



'tip I 27ciknn / 



^{k + 2)(Mk + 2) + 2 



To summarize, the Z/e-parafermions coupled with one J7(l), with symmetry generators Jo ~ idcj), 

J\+M/2 ~ , 0e^ M+fe( ^, and Ji +M / 2 ~ ^ e V M+ fe^ 5 will have ( fe + 1 )(^ ffe + 2 ) multiplets. We can use 
the (£; n) to label them, with £ the parafermion loop the state sits and together with n determines 
the U(l) charge assignment. When k is odd (even), £ = 0, . . . , (£ = 0, ...,| + 1), n = 
0, . . . , Mk + l(n = 0, . . . , ^ if £ = |). 
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